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Abstract. 

The sigma form of the Painleve VI equation contains four arbitrary 
parameters and generically the solutions can be said to be genuinely "nonlinear" 
because they do not satisfy linear differential equations of finite order. However, 
when there are certain restrictions on the four parameters there exist one 
parameter families of solutions which do satisfy (Fuchsian) differential equations of 
finite order. We here study this phenomena of Fuchsian solutions to the Painleve 
equation with a focus on the particular PVI equation which is satisfied by the 
diagonal correlation function C{N,N) of the Ising model. We obtain Fuchsian 
equations of order N + 1 for C{N, N) and show that the equation for C{N, N) 
is equivalent to the A^*'' symmetric power of the equation for the elliptic integral 
E. We show that these Fuchsian equations correspond to rational algebraic 
curves with an additional Riccati structure and we show that the Malmquist 
Hamiltonian p,q variables are rational functions in complete elliptic integrals. 
Fuchsian equations for off diagonal correlations C{N, M) are given which extend 
our considerations to discrete generalizations of Painleve. 
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1. Introduction 

The correlation functions of the Ising model were first calculated by Kaufman and 
Onsager ^ in terms of determinants whose elements are certain hypergeometric 
functions. For this reason it follows from a theorem on holonomic functions j2| that 
they must satisfy linear ordinary differential equations. However, these correlations 
also have a remarkable connection with nonlinear equations as well. The first such 
result was the expression as T — > of the scaled correlation function in terms of a 
PHI function by Wu, McCoy, Tracy and Barouch j3j in 1976. Subsequently in 1980 it 
was shown for arbitrary fixed T by Jimbo and Miwa @] that the diagonal correlation 
C {N, N) is given in terms of a P VI function and by McCoy, Wu ^ and Perk |S] that 
the correlation at a general position C{M,N) and its "dual" C*{M,N) satisfy some 
remarkable quadratic identities, or double recursions which are discrete generalizations 
of the Painleve ODE's. 

The Painleve representation of the correlation functions is by now well known but, 
curiously enough, almost nothing is known about the corresponding linear equations 
beyond the fact that the diagonal correlation function C(l, 1) is a particular case of 
the hypergeometric function. In this paper we will study these linear equations for 
the Ising correlation functions and the much more general question of when solutions 
of the PVI equation will satisfy Fuchsian differential equations. 

The most general four parameter dependent sigma form of Painleve VI can be 
written as [71 E] 

C{t{t - 1)C")' + {2C{tC - C) - C" - V^V2V3V4? 

= (C + vDiC + vl) iC + vl) iC + vl) with: (1) 

C = t(t-1)^^^^^ +Kit +K2 where: (2) 
dt 

Ki = - - and: (3) 

K2 = — -(wiW2 — — i'iU4 — ^2^3 — ■y2i'4 + (4) 

This is a second order nonlinear equation which allows branchpoints only at the three 
points t = 0,1, 00 and locally near these singularities the function t has, following 
Jimbo's expansions ^j, an expansion of the form 

00 00 

T = xP^ x'^'+'^^Y'i]-'^ ■aj{n,k;a)-x'' (5) 

k—~oo n—Q 

where x is the local variable at t = 0,1, 00, and two boundary conditions for the 
second order PVI equation specified by a and 77 will in general be different at the 
three singularities. The coefficients aj{n, k; a) depend on the value of j = 0, 1, 00 and 



satisfy aj{n, — fc, a) — aj{n, k, ~a) and we note that 

Pa ^ {a"^ ~ {vi + V2 - V3 - Vi)^} / A (6) 

Pi= {a^-{vi+V2-V3 + Vi)^}/A (7) 

Poc = ay A + Ki (8) 



Comparison of Q with the well known expansion of Jimbo IHj reveals that many of the 
coefficients in Jimbo's expansion vanish identically. Several aj{n,k,a) are explicitly 
given in Sec. 2.1. 
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In general the local expansion |(SJ) has an infinite number of confluent singularities 
which indicates that it cannot satisfy a linear differential equation. Therefore even 
though the most general solution of the PVI equation cannot satisfy a linear equation, 
the specific boundary conditions which specify the solution to be the physical diagonal 
correlation function of the Ising model will allow a Fuchsian equation of order 
generically greater than two to be satisfied. 

In this paper we study this phenomena of the existence of boundary conditions 
for which solutions of certain PVI equations satisfy Fuchsian differential equations§. 
There are several ways in which this phenomenon may occur. One way is that 
conditions can be found on the four parameters Vk and on a such that the general local 
expansions at t = 0, 1,cxd degenerate by having the coefficients aj{n,k;a) all vanish 
if k is sufRciently large. This will give a one parameter family of solutions which has 
only a finite number of confluent singularities. We study this mechanism in detail 
in Sec. 2.1. However, there may also exist one parameter families which cannot be 
obtained from the two parameter families lO by specialization. An example of this is 
given in Sec. 2.4. 

For concreteness we will consider in detail the specific PVI equation for the 
diagonal Ising correlation obtained by Jimbo and Miwa 0]: 

t{t^l)<7"Y = (9) 



((i- l)o-' - aY - 4cr'((t- l)a' - a- 1/4^ (ta - cr) 

which is obtained from by setting 

= = N/2, t-2 = (1 - N)/2, V3 = {1 + N)/2 (10) 
a = C +Nh/4:-l/8 (11) 
The diagonal Cn = C{N, N) is related to ct for T > by 

ait) = <(t-l). Aiog(C^) -1/4 

with t^ (smh{2Jy / kT) ■ smh{2 Jh/kT)y < 1 (12) 
where Cat ~ t^/^ ^^en t 0, and for T < by 
ait) = t(t-l)-^log(Cjv) -t/4 

with t = fsinh(2J„/fcr) • sinh(2 J,,/A:T)) ^ < 1 (13) 



where Cn — 1 when t — > 0, and where the variable Jy iJh) is the Ising model vertical 
(horizontal) coupling constant. The detailed specification of the behavior of ait) near 
t — Q needed to uniquely specify Cm as the diagonal Ising correlation function are 
sketched in section ()2.4|l (see for instance equation (|53|l '). Do note that since all the 
calculations of this paper are systematically checked with high-temperature expansions 
when available, we introduce a variable t which is the inverse of the one of Jimbo 
and Miwa For integer N the equation © is in the class of so called "classical" 
equations 8^ which arc known to generate Toeplitz determinants whose elements are 

§ For a warm-up on Painleve VI, sigma form of Painleve VI, and on the question of the holonomic 
solutions inside Painleve VI we recommend two magnificent papers in French, one by Garnier 
and the other one by Okamoto 21 {^^e also in English llll 'l. 
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hypergeometric functions |7||Sm]. We present the Fuchsian equations satisfied by 
Cn for small values of N in Sec. 2. These equations have remarkable structure 
and in Sec. 2.2 we show that the associated A'' + 1 order differential operators are 
homomorphic to the 7V-th symmetric power of the second order differential operator 
associated with the elliptic integral E. In Sec. 3 we present an algebraic formulation of 
the Fuchsian equations for C{N, N) by studying the Riccati formulation of solutions to 
PVI for iV = 1, 2 which are related to differential structures on certain rational curves. 
In Sec. 4 we extend our considerations to the discrete generalization of Painleve 
VI, namely a quadratic double recursion on the two-point correlation functions 
C{N, M) together with their dual C*{N, M). We will show that these structures 
can be generalized, mutatis mutandis, to the C{N, M)'s. The C{N, M)'s are also 
solutions of Fuchsian linear ODE's, with a quadratic increasing order. The associated 
differential operators are now homomorphic to direct sums of iV-th symmetric power of 
the second order differential operator associated with the complete elliptic integral E. 
The C(iV, M)'s are actually sums of several homogeneous polynomials in the complete 
elliptic integrals E and K. This is a consequence of various remarkable simplifications 
in the "discrete Painleve" double recursions, like the fact that algebraic or rational 
expressions become polynomials by remarkable factorizations and by the occurrence of 
perfect squares. Combining these various results together, one has some quite curious 
and fascinating alchemical wedding between complete elliptic integrals, rational curves 
and discrete generalizations of Painleve VI (and Hirota-Backlund transformations). 
The confrontation between the non-linear Painleve world and the linear Fuchsian 
world (Painleve versus Fuchs) yields the emergence of quite interesting structures of 
differential nature but also of algebraic geometry nature. We finally see in Sec. 5 that, 
in the case of the C(7V, N) holonomic solutions, the p and q Malmquist's variables 
corresponding to the Hamiltonian structure of the sigma form of Painleve VI are 
remarkably rational expressions of E and K , and even rational expressions of E/K . 
We have the same result for the a and a' variables. These last results are in complete 
agreement with the previous mentioned results, namely the rational character of the 
algebraic curves corresponding to the existence of holonomic solutions C{N, 7V)'s for 
the sigma form of Painleve VI, and the existence of simple Riccati equations for the 
uniformizing parameter. 

The number of new exact results we have obtained being quite large and the 
explicit formulas for some of these results being quite cumbersome, we will just sketch 
here these new exact results, giving the simplest formulas. More exhaustive formulas 
will be given in forthcoming publications. 

2. Solutions of sigma form of Painleve VI and Fuchsian linear ODE's 

We consider, from now on, the isotropic square Ising model and the high temperature 
regime, i.e., t — s'^ where s ~ sinh(2J/fcr). The introduction of these two variables, 
t and s, may look a bit redundant: the variable t is well-suited to write down our 
results on diagonal correlations functions, while the variable s is clearly better suited 
for non diagonal correlations. The results for the low temperature regime are similar. 
The diagonal two-point correlation functions of the square Ising model C{N, N) and 
its dual C*{N,N) can be calculated from Toeplitz determinants [Tl 1121 IT^ : 





(14) 
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C*(Ar,iV) = (-l)^dct(a,_,_i), l<t,j<N (15) 
where the a„'s read in terms of 2F1 hypergeometric function 

" - ^"^/^■^t"/2+i/2.2^^,(i/2,„+l/2;n + 2;t), n>-l(16) 



(n+ 1) 
(l/2)-„-i 



gi'"/'-^/'-2Fi(-l/2,-n-l/2;-n;t), 



n < -1 



- (-n-1)! 
where (a)„ is the usual Pochhammcr symbol. 

The diagonal two-point correlation functions of the square Ising model C{N, N) 
and C* {N, N) being given by the Toeplitz determinant p4|l whose entries are solution 
of linear second order differential equations, they are necessarily solutions of a linear 
differential equation, with order A'^! • 2^ as an upper bound for generic entries of the 
determinant. 

Since the diagonal two-point correlation functions of the square Ising model 
C{N,N) are given by the determinants H14|l . it is straightforward to obtain a 
sufficiently large number of series coefficients and to get the linear differential equations 
satisfied by these series f^ElEl- Denoting by Dt the derivative with respect to the 
variable t, the first linear differential operators Lat^v corresponding to the C{N,N) 
are 



n (t-2) ^ 1 1 t + 2 



4 i-5 ,1 (41- llt-2<2) 
L33 = Df +2 -j^ ^ ■ Df +- ^ ^ — • D? 

* it-l)t * 2 t^t~lf * 

1 (2t^ + 2t-5) 9 15+13t + 4t'^ , , 

20 , (lU+7t-2t^) 3 
L44 = • Dt^ + ^ ^ ■ Dt^ 

t{t~ 1) it - if i2 



1 (322 + 95i-9i2 _ 16*3) 

2 (i-l)'<3 



ZT. Dt^ (20) 



(97 40 1- 10*2 -12*3) 32-t-33t-h20*2 + 5t3 
+ -^^ ^ ■ Dt - 4 T , 

[t-lff^ (i- 1)^*5 

(*-h7) ,1 (52*2 + 483*+ 1617) 



tit -I) 4 (*- 1)2*2 

• Dt^ (21) 



1 (4*3 + 370*2 + 1707* + 3503) 

2 (-l + *)^*3 



1 (1552 *4 - 1016 *3 - 13191 *2 - 29618 * - 29855) 
16 *4 (*- 1)* 

5 (720 ** + 640 *3 - 2175 *2 - 6912 * - 88OI) 



Dt 



2 



16 (*- 1)^*5 



Dt 



25 784 *4 + 3428 *3 + 6921 *2 + 8650 * + 7865 
+ ^ ^ -5 



64 (* - if *6 
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(4 + <) « (81 + 39i + 7i2) 



where 



--^^•A-^^V, (22) 

Ni = 10162 + 7059 1 +241U2 + 376 1^, 

iVa = 37973 + 35162 1 + 17893 + 5116 + 500 

iVs = -28706 - 55327 46180 1^- 21437 1^- 3358 + 1736 ^^ 

iVi = -390548- 402496 <- 240997 i^- 63239 

+ 24152t'' + 25088 <^ (23) 
iVo = 23814 + 26839 1 + 24583 f + 16599 + 7345 + 1620 

These operators are of order + 1 and are irreducible. We further note that, 
in contrast to the Fuchsian equations for the n-particle contributions ^("^'s of the 
susceptibihty of the Ising model [21 ^1 1161 ITTj , the Fuchsian differential equations 
satisfied by the C(iV, A^)'s have no apparent singularities. The linear differential 
operators, L^^^, for the C* {N, N)'s are obtained by the change t into 1/t in the 
previous differential operators. 

These Fuchsian differential equations H17I22|1 for the C {N, N) 's have the following 
general form : 



N+l 

E 

i=2 



t\t~iy-^P^^''\t)-D\+t{t-l)Pi'^\t)-Dt+P^'^\t) ^ (24) 



where P/^' [t) is a polynomial in t of degree A^+1 — ifori = 2, •••,A''+1 and P['^^ (t) 
and PQ^\t) are of degree iV — 1. 

The only singular points of H24(l are the three regular singular points t — 0, oo, 1. 
From the indicial equation of the differential equations for the first Ljvat's, we infer 
the remarkably simple expressions of the critical exponents p^°°'> and p^^'^ at 

respectively the regular singular points t — 1, t = oo and t = 

= (n-iy (25) 
p(°) = -i + ^iV+i(n+l)(n + 2) (26) 



1 (-1)" f-D" 



where n = 1, 2, • • • + 1. 



2.1. Local solutions at t = 0,1, oo 

It is of interest to compare the local expansion ||SJ) of the PVI equation with the 
exponents of the Fuchsian equations! |. For concreteness we concentrate on i = 1~ 

II Recall that, for Ising case and for T > Tc, t = t^l^ Cn (resp. r = Cat for T < Tc). 
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which corresponds to T = r+ in the Ising model. We have the foUowing coefficients 
in © valid for < a < 1 

ai(0,0;a) = l, (27) 
ai(0, —1; a) ~ ai(0, 1; —a) 

= -.^ ^ {a ~ vi - V2 " V3 + V4){a - vi - V2 + - Vi) x 

loa^(l — a)^ 

{a — vi + V2 — V3 — Vi){a + vi ~ V2 ~ — V4), (28) 
1 

ai(l,0;a) = — — + -{-V1V2 + uiWa + V1V4 + V2V3 + V2V4 - V3V4) 
o Z 

+ TT^ (^1 +V2+V3- V4){vi +V2-V3+ V4) X 
8q^ 

(wi - U2 + W3 + W4)(wi - W2 - W3 - ^^4), (29) 
a.(0,-2;a) ^ a.(0,2;-a) ^ ,^5(0 -"iWa- 2^ (, - 3)^ 

X [(a - 2)2 - (wi + V2 +V3- V4f] [{a - 2f - {vi + V2 - V3 + V4f] 
X [{a - 2)2 - {vi -V2+V3+ V4f] [{a - 2f - {-vi +V2+ V3 + Vif] 
For the Ising case (|1U|I this reduces to 

Pi = aV4, (31) 

ai(0, -1; a) = ai(0, 1; -a) - ai(l, 0; a) = ~ " \ (32) 

loa 8 

ai(0,-2;a) = oi(0,2;-a) = 256 (a - 2f ^ ' 

When used in ||SJ) these expressions will reproduce the N +1 exponents of L mn at 
t = 1 where, in the limit a ^ 0, the terms in (O with become {t—1)^ ln'^(<— 1). 

We see from that, when a = 0, ai(n, ±2;0) = for iV = 1 which is consistent 
with the fact that C(l, 1) satisfies a second order linear differential equation. We have 
carried the expansion to order {t—lY"^. In particular we have obtained the coefficient 
of {t — 1)^ In'^ 1^ ~ 1| and have verified that it vanishes for = 1, 2 and have obtained 
all terms in the expansion of C{N, N) given in |T8j . 

More generally the conditions that there exists a value of a such that aj(n, fc; a) = 
for all k sufficiently large is a condition necessary for r function of the PVI equation 
to satisfy a linear differential equation of finite order and the series 



oc 



aj{n, k; a) ■ x'^'' (34) 

n=0 

will be solutions to the Fuchsian equation. For example one condition for a second 

order Fuchsian equation is ai(0, l;a) — 0, ai(0, — 2;a) = 0, which are satisfied if, 
respectively, 

a = —vi — V2 + V3 — V4, a — 2 = —vi + V2 — V3 — V4 (35) 

implying V2 ~ V3 = —1, which is the restriction Forrester and Witte B needed for 
a solution of PVI to satisfy a hypergeometric equation. This condition implies that 
the r functions are determinants of hypergeometric functions. We thus see that, at 
order a;P+4+(i-a)^ the local expansion provides a necessary condition for the reduction 
of a one parameter family of solutions to PVI to a solution of a second order linear 
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differential equation. By examining the vanishing of a(0,k]a) for higher values of 
k necessary conditions for the existence of one parameter families satisfying higher 
order linear differential equations will be obtained. Similar necessary conditions can 
be obtained from the local expansions at t = 0, oo. 



2.2. The Fuchsian differential operators as N-th symmetric power 

The most profound and surprising structure of the solutions of PVI which satisfy 
Fuchsian equations is, however, not seen in these local expansions and, thus, it is 
important to observe that the operators L^^n given in (|17I22|) for C{N,N) have the 
remarkable property that they are equivalent [I to the A^-th symmetric powerff of 

An ■ Lnn = Sym^{Lii) ■ Rn (36) 
The first and Rm intertwinncrs read for N = 2 : 

2^2 l(31i-23)< ^ 3 15i-7 

A2 = t^Df + - ^ '— -Dt -\ 37 

*4i-l 4 t-1 ^ ' 

t'-D^+lt-D, (38) 

We have calculated exactly these intcrtwinners up to iV = 6 but the expressions are 
too large to be given here. As a consequence of this property (|36|l the differential Galois 
group oi Lnn is not a SL{N + 1, C) group as we could expect at first sight, but an 
SL(2,C) group in the symmetric power representation. We expect that this property 
extends much more generally to other solutions of the general four parameters PVI 
which satisfy Fuchsian equations. 

Another consequence of 1)36(1 is that the solutions of this order N + 1 differential 
operator Ln n are actually homogeneous polynomials of degree N in the two solutions 
of ill, see IlIIEllESl. 

Let us now introduce the two elliptic integrals 

K = 2F1 (1/2, 1/2; 1; s^) , E ^ 2F1 (l/2, -1/2; 1; s^) (39) 

and the second order linear differential operator for E (Ds denotes the derivative with 
respect to s): 

Le - +— -4^^ (40) 
s s^ ~ I 

This operator actually identifies with LJ^- 

One can easily show that the second order linear differential operator Ln 

(associated with C(l, 1) and written in the variable s) and the second order linear 

differential operator Le are equivalent : 

• As + 6s2) . ill = Le ■ • Ds - 2/s') (41) 

More generally one can show in the s variable, that the Lnn^s are actually equivalent 
to the L*j^i^^ s. Since K can be simply expressed in terms of E and its first derivative, 
the Cn.n's are thus solutions of an operator which is homomorphic to Sym^ {Le) ■ 

An ■ Lnn — Sym^iLE) ■ Rn, or: (42) 

Lnn -Bn = Sn ■ Sym^{LE) (43) 

(t For the equivalence of differential operators see (e.g.) I19ll2(jllljll . 

ft For the definition of the symmetric power of a differential operators see (e.g.) 1211 1221 |23I . 
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where the intertwinners -Bjv and Sj\j (or Ajv and Rn) are Unear differential operators 
of order iV. In fact, beyond C{N,N), relations (jSHl, ESI relate all solutions 

of Lnn to Sym^ {L e) ■ From (|43|) one can easily deduce that the diagonal two-point 
correlation functions C{N^ N) can be deduced as the action of a linear differential 
operator of order N on the A^-th power of the complete elliptic E : 

C{N, N) = Bn{E^) (44) 

The expressions of the intertwinners Sat can also be retrieved from the 
determinental expressions of the C(N,N) and the relations between the hypergeometric 
functions a„, a„_i (see ((T^ . (fTBI) ) and its derivative. This gives a general method to 
obtain the differential operators L^ js[. 

2. 3. The C {N, N) 's as homogeneous polynomial of the complete elliptic integrals E 
and K 

The property (|36|l . or H43|l can be illustrated by considering the specific solution 
C(N, N) of the iV + 1 order differential equations Lnn- The matrix elements a„ of 
the Toeplitz determinant representation may all be expressed as linear combinations 
of the elliptic integrals E and K, and, thus, C{N,N) will be given as polynomials 
in these functions and this is in agreement with the previous relation H44|) . For low 
orders these polynomials have been presented by Ghosh and Shrock 24 . For example 

C{2, 2) - ^ • (3 {s^ - if + 8 {s^ -l)-EK - (s* _ 5) • E^) 
C{3, 3) = ■ P^iE, K), where: P^iE, K) - 

(335^-1) {s^-lf-K^ +3 (s8 + 48/-l) {s^-lf-EK^ 
-3 {s^^ + is^-ms^ + l)-E^K (45) 

- (1 + 21 s^ - 96 s'^ + 10 5^2) . E^ 

We note that these expressions are respectively quadratic and cubic homogeneous 
polynomial in E and K. We have obtained similar expressions for all the C{N, N) 
and C*{N, N) for N = 4, 5, 6, • • • , 21, and relation (|44|l gives similar relations for 
any values of N. They are homogeneous polynomial of degree N in the complete 
elliptic integrals^ E and K, with simple rational coefficients (a polynomial in s 
with integer coefficients divided by some power of s). From a physics viewpoint one 
should note that the particular rational coefficients one gets in front of the monomials 
. K^^^ , are far from being arbitrary as a general formula like H44|l could suggest. 
These coefficients are such that, for instance, the linear differential equation for the 
C{N, N)^s has no apparent singularities. Furthermore, the contribution associated to 
the various monomials E'^ ■ K^~^ clearly have poles {s~^^ or in the previous 
example (|45|l ). These coefficients are also "fined-tuned" in such a way that, for 
instance, these various poles cancel together, in order to give an expression with a 
well-defined high-temperature series expansion (series at s = 0). We have many 
other remarkable properties corresponding to the behavior of the C{N, N)^s near 
s — 1 or s = 00. 

^ This result can also be found in the Eqs. (2.16)-{2.19) of 1251 . which also show very explicitly that 
C(N, N) is a. homogeneous polynomial of E and K of degree A'^ for all A'^ (something that is already, 
albeit less explicit, in the appendix of MontroU, Potts, and Ward 1121 .'). 



Painleve versus Fuchs 



10 



2.4- Non-trivial disentangling of solutions of linear Fuchsian ODE's near t = 0. 

Let us make here a comment on the existence of surprisingly simple hypergeometric 
solutions of the A^-dependent sigma form © of Painleve VI. Consider the second 
order differential operator: 

which has regular singularities at t = 0, t = 1 and t ^ oo with respectively the critical 
exponents (±iV/2), (1/4,1/4) and (l/4± Af/2). 

It can be verified that any linear combination of the two solutions of (|45|l satisfies 
the A'^-dependent sigma form (j^ Painleve VI equation /or arbitrary N, not necessarily 
an integer. For instance, when N is not an integer, one has the two following solutions 
of ®: 

a = t{t-l) - - where : r = /+ + A • /_ (47) 

where f± are the two independent solutions of (|46|l : 

/± = t±^/2.(l-t)i/4.2Fi([l/2,l/2± iV], [1 ± iV],t) (48) 

When the parameter iV is an integer (and only in this case), that is to say in the 
Ising case we are interested in, the second order differential operator Lh is, after 
conjugation by (1 + s^)^/^, equivalent to Le', when N is an integer, one solution is 
given above in term of a hypergeometric function analytic at i = 0, and the other one 
has a logarithmic singularity aX t — Q (and similarly for t = \ and t = oo). 

At first sight the existence of such "additional" solutions should not be seen as a 
surprise: we certainly expect the solutions of the A'^-dependent sigma form of Painleve 
VI that are also, at the same time, solutions of a linear (Fuchsian) ODE, to be a quite 
complicated "stratified" space. However, let us focus on the series expansion a.t t — 
of the analytic solution of H4t)ll , which simply reads 

_ 1 r(2jv + i) 
^ " 4^r(7v + i)2 ■■^+ 

= co(A^)-i^/2^ci(iV).t^/2+i+c2(iV).t^/2+2 + ... (49) 
The coefficients Ck{N) in the series expansion of H49|l read : 

CkiN) = (50) 

■ 3F2 ([1/2, 1/2 + N, -fc], [1 + N, 5/4 - k], 1) 



1 r(2A^+l) (-1/4), 



4^r(iV + l)2 fc! 

Let us now consider the series expansion of the diagonal correlation functions 
CiN,N): 

C{N,N) = do{N) ■ t^'"^ + di{N) ■ + ^2(iV) • + • • • (51) 

where do{N),di{N) and d2{N) read respectively : 

r(2iv + i) j_ r(2jv + 1) TV 

r(iv + i)r(Ar + 1) 4^' r(7v + i)r(7v + 2) 4^+1' "' ^ ' 

One has the following result, that may look quite surprising at first sight: the 
coefficients Ck{N) of the solution and the coefficients dk{N) of the diagonal 
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two-point correlation functions C{N,N), solution of the order + 1 Fuchsian ODE 
are identical up to k — 3N/2 + 1: 



The coefficient in (|53f) in front of ^3^/2+2 (^aj^ i-,g gggj^ a,s the initial condition defining§ 
C{N, TV). Seeking for conditions allowing solutions of the sigma form of Painleve VI to 
be also (the log-derivative of) solutions of linear Fuchsian differential equations, this 
difficulty to disentangle, near i = 0, a solution of a second order differential equation 
and a solution of linear Fuchsian differential equations of arbitrary A'^-l- 1 order, seems 
to indicate that series analysis like © may not be the easiest approach to take into 
account such subtlett fine-tuning: we need a less analytical and more "global" algebraic 
approach. 

3. Algebraic viewpoint of the Fuchsian differential equations 

The existence of C{N, N) as solutions common to the sigma form of Painleve VI 
equation and to linear Fuchsian differential equations can be addressed on an effective 
algebraic geometry approach of differential equations as introduced explicitly by J.F. 
Ritt |27ll28j . This approach amounts, when working with various linear and non- linear 
differential equations, to introducing as many variables as the number of derivatives 
of the function we study. The analysis of the compatibility between these various 
linear and non-linear differential equations will correspond to considering an algebraic 
variety given by various polynomial relations on these variables. These relations can 
be studied from the algebraic viewpoint (parametrization when the genus is zero or 
one, birational transformations^, singularity analysis, blow-up, etc.). The very last 
step, recalling that the various introduced variables are not independent but can be 
deduced from each other by successive derivation, provides further constraints. In 
other words a set of differential equations is seen as an algebraic variety plus some 
differential structure on top of it. 

Let us show how this algebraic viewpoint of differential equations works in our 
(subtle) compatibility problem of the sigma form of Painleve VI and the Fuchsian 
linear ODE's of arbitrary order N + 1. The correlation function C(l, 1) satisfies 
a second order linear differential equation which can be written in a Riccati form in 
terms of a{t) and cr'(t). More generally, the N+1 order Fuchsian linear ODE satisfied 
by the C{N, N)'s can be written in a "generalized Riccati form [23 ED]" in terms of 
a(i), a'{t) and its successive derivatives a^^\t) up to n = N (where a{t) is deduced 
from C{N, N) by the logarithmic derivative relation (|12|l ). Similarly, the sigma form 
of Painleve VI equation is not seen as a non-linear ODE, but as a polynomial 
relation between the three variables cr(t), cr'(t) and <j"{t). 

Introducing the variables 5*0 = a{t), Si — cr'{t), S2 = f"(t), etc., the third 
order Fuchsian linear ODE for C(2, 2), yields a "generalized Riccati form" which is a 

§ For T < Tc, one has C(Af, A^) - (1 - i)!/^ = 1/4 . {(1/2)^ (3/2)iv)/((Af + 1)!)2 ■ t^+^ +■■■, 
the coefficient in front of t^'^^ corresponding to the initial condition defining C'{N, N) in the low- 
temperature regime 1261 . 

(t Cauchy's theorem does not apply to PVI at t=0 or t=l. As a consequence, even with given 
boundary conditions (a large set of first terms in the series), there can be "branching" in the series 
computation. These subtle "branching" series calculations will be adressed elsewhere. 
^ At this step it is worth recalling that Biicklund transformations are actually birational 
transformations in "some" variables. 




(53) 
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polynomial relation between 5*0, and ^2 

64 {t-lf S2-l6t{8t + 5) {t-l)-Si 

+ l92t{t-l)SaSi+6ASo^ - 16 (16i+ l)S'o^ 

+ 4 (32<2 + i6i-2l) 5o + 45 = (54) 

The elimination of the variable 5*2 between this "generalized Riccati form" and 
seen as a polynomial relation between the three variables 5*0, Si and 5*2 yields an 
algebraic relation between 5*0 = a{t) and Si — cr'{t) which reads: 

(4S'o - 3) (6450^ - 16 (16 1 + 1) + 4 {Mt^ - 16t- 21) • S'o + 45) 

- 32 1 (4 5*0 -3) (t - 1)(8< - l-4S'o) •S'l 

+ 256t2(i- l)^^!^ = (55) 

which is compatible with (|54ll and (O . This can be checked by eliminating S2 between 
the derivative of H55|) and H54I) or to get again (|55|l . Or directly by plugging a series 
expansion or an exact expression of C(2, 2) in (|55|l . 

Seen as a relation between Sq and Si (the variable t is considered as a simple 
parameter), the algebraic curve (|55|l is actually a rational curve. It can thus be 
parametrized in term of two rational functions: 

3 A2 • + Ai • u + Ao 

4 B2 ■ + Bi ■ u + Bq 



Si 



3 {ai-u + ao) ■ {C3 ■ + C2 ■ u"^ + Ci ■ u + Co) 



where: 



t {B2-u^ +Bi-u + Bof 

ai = -6 i - 3 + 8 ao = 4 • (1 - 2 1) 



Aa = -176 + 48^- 320<2 + 256^^ 

Ai ^ 120+ 184i-144i^ + 768t^-512i*, (57) 

A2 = 9-57^ + 24^2 + 76<^-448^'' + 256t^ 

Bo = 192 f - 272t- 112, Bi = -8 (3i + 1) {l6t^ - 26i- 3) , 

B2 = 45 + 51 < - 168 - 260 + 192 1^, 

Co = 1088 + 384 i + 2624 + 1280 1^- 1536 1^, 

Ci = -1296 - 2816 1 + 688 i^- 7776 1^- 3840 + 4608 ^^ 

C2 = 108 + 1848 1 + 636 f - 3328 + 8304 i'' + 4416 - 4608 1^, 

C-i ^ +189 + 36 1 - 1323 f + 210 f"^ + 2460 - 2792 



- 1856i^ + 1536i^ 

In the spirit of the "algebraic viewpoint of differential equations" [23 EH| , having 
performed the algebraic geometry calculations we had in mind, we now recall that there 
is some differential structure on this rational curve by imposing that the variable Si 
is actually the derivative with respect to t of the variable So'- 

_ dSo _ dSo du dSo 
^'-^-^'A^^ ^^^^ 

yielding, after some quite nice simplifications, that ^ is not a rational expression of 
u, as one could expect at first sight, but a quadratic polynomial in u, which gives a 
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simple Riccati form: 

16i(t-l)(6t2-5i-9) ~ = 

(63 - 135 i - 120 - 140 + 192 1^) ■ (59) 
+ 8 (15 + 51 < + 46 - 60 1^) ■ u - 272 - 112 i + 192 1'^ 

that can easily be associated with a linear second order differential equation bearing 
on some function F: 

_ 1 dF _ 1 192 _ 140 13_ 120^2 -135^ + 63 
" ~ F'~dt ^ "16 t{~l + t) (6i2-5t-9) " ^ ^ 

Similar calculations can be performed for = 3, the generalized Riccati form 
for the Fuchsian linear ODE of order four is now a polynomial relation of the form: 

^3 = P{So, Si, S2; t) (61) 

where P is a polynomial of the three variables Sq, Si and S2, the coefficients being 
rational function (with integer coefficients) in the variable t seen as a parameter. In 
order to combine§ this generalized Riccati form H61() with lO for N —3, we need, in 
order to perform eliminations of variables (ideal of polynomials), to rewrite 0, the 
sigma form of Painleve VI taken for N = 3 as a relation between a, u' and a" and 
a^^^ as well. This is easily obtained by performing the derivative of lO with respect 
to t, thus getting a polynomial relation between cr, a' and a" and cr'^^. Considering 
this last polynomial relation and the generalized Riccati form (|61|l . we can easily 
eliminate 6*3 — a^'^^ , getting a new polynomial relation on iSq, S\ and 15*2. We can, 
now, eliminate 6*2 between this new polynomial relation and @ for = 3 which is 
also a polynomial relation on 6*0, S'l and 5*2, in order to get, finally, a polynomial 
relation on 5*0 = a and S'l = a' only. This final relation reads: 

4096 [t - \f ■ Sf + 256 {t - if Q2 ■ S^ (62) 

-16t(t-l) gi-Si - (45-8 (2i + 7) So + ieSo^) -Qo = 

where: 

Q2 = 48S'o^ -8 (22t+13)-S'o + 55 +448t +64^^ 

Qi = -768S'o'^ + 256 (22t + 13) - 32 (376^^ + 584i + 125) Sq^ 

+ 16 (384 + 1984 + 766t + 25) Sq + 1 125 + 2880 1 ~ 25920 
Qa = 1575 + 16 (576 _ no t - 145 - 96 1^) So (63) 

- 32 (56 1 - 9 + 264 1^) Sq^ + 256 (10 1 + 3) Sq^ - 256 Sq'^ 

Similar calculations (of ideal of differential equations seen as ideal of polynomials), 
can be performed, mutatis mutandis, for = 4, 5 and 6. These eliminations yield 
polynomial relations in t, Sq — a and Si = a' of the form: 

i=N 

J2 f{t-lTn{So,t)-Si' = (64) 

i=0 

where the Pi{So, tys are polynomials in t and 5*0 = cr, of degree 2i in Sq- Again, 
these relations (|64|) seen as algebraic curves in Sq and Si {t being seen as a parameter), 
are rational curves. From the previous remark that the C{N, N) are homogeneous 

§ Or, in mathematical wording, to calculate the ideal of these two differential equations. 



Painleve versus Fuchs 



14 



polynomials of E and K one can easily deduce that Sq — a and Si = a' are rational 
expressions of the ratio r — E/K (or E' /E). 

Now, similarly to the previous calculations, recalling that the variable 5*1 is the 
derivative with respect to t of the variable S'o, one also finds Riccati equations similar 
to (|59() for the uniformizing parameter u: 

$ = P2{t)-u'+ Pi{t) ■ u + P^{t) (65) 
at 

where (3o{t), Pi{t) and /32(t) are quite simple rational expressions of t, the Riccati 
equation (|65|l having only i = 0, t = 1 and t = oo as regular singularities. The 
calculations are too large to be given here and will be detailed in a forthcoming 
publication. 

Note that, in such "global" Riccati algebraic approach, one has to be careful 
because of the existence of many singular || solutions of corresponding to algebraic 
T functions : 

a - t(t-l)-^log(T) -1/4 (66) 

r = e -{l-tf, [AP-lf N'^ + lQl3{'ia + l){a + P) = 

like, for instance, (a, /?) being {-N/2, -1/4 • N/{N - 1)), (-1/8 ■ {AN^ + 1), N'^) 
or (—1/4, 1/4), and especially {N/2, — 1/4 • A^/(iV+ 1)) which corresponds to a series 
expansion with leading order similar to H49|l . 



4. Generalization to non-diagonal correlation functions C{N, M) 

Most of the results, previously displayed, can be generalized to the non-diagonal 
correlation functions C{N, M) of the square Ising model. The C{N, Af)'s are also 
given by determinants (see \12\) whose entries are holonomic quantities solutions 
of linear differential equations of order three. The C{N, M)'s are thus holonomic 
solutions of linear differential equations. At first sight the growth of the order of the 
corresponding differential operators should also be exponential in TV and M. 

We found that the order of these linear differential operators is, again, not growing 
exponentially with N and M but has a quadratic growth order and depends on the 
parity of M — N. For all the Fuchsian linear differential operators we have obtained 
{N and M < 6), the order can be reproduced by : 

q= i-(M + Ar + 2)-(4 + (3- (-l)*^-^)-|A/-iV|) (67) 

These linear differential operators L^m are too large to be given explicitly here. Let 
us just give one of them, namely the linear differential operator L12, corresponding 
to the simplest non-diagonal (and non horizontal or vertical like C(0, N) or C{N, 0)) 
two-point correlation function. The linear differential operator L12 reads 

5{2s^ + 3)Dt , q,-D^, 



L12 - 

+ - 



S(l + S2) s'^{l + sf{l- Sf{l + S-^Y 

92 • Dl , qi-D, 



(l-hs)3(l-s)3(l + s2)3 s\l + SY{1- SY{1 + S-^Y 



S5(l + S)3(l-S)3(l + S2)5 
II We use here the terminology of singular solutions of differential equations 1271 1281 . 



(68) 
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where the polynomials Qi read : 

53 = 13 + 30 s^- 78 s''- 50 + 53 (69) 
(72 = 5s^2 - 7si° + 345^ - 128s^- 655-*- 97 5^ + 2 
qi= - 5 s^-* + 2 _ 67 fiio - 118 - 816 + 157 - 76 - loi 
qo= - 192 s^" + 1840 s^ - 453 + 127 - 15 - 27 

Let us comment on the remarkable simplifications we encountered when com- 
puting the C(-/V, M)'s from the quadratic double recursions (discrete generalizations 
of Painleve equations) they satisfy ^Hl together with the C*{N, Af)'s . From the 
expressions of the C{N, N)'s as homogeneous polynomial in E and K, and the ex- 
pressions of C(0, 1), we can obtain the C{N, M) and C*{N, M), step by step using 
this quadratic double recursion At first sight these C{N, M)'s should be given 
as rational expressions of E and K and, in some cases, as roots of quadratic poly- 
nomials with polynomial expressions in E and K. Remarkably, as a consequence of 
factorizations and simplifications in the numerator and denominator of these rational 
expressions, and the occurrence of a perfect square in the case of roots of quadratic 
polynomials, the C{N, Af)'s are actually always given by polynomial expressions in 
E and K , that are no longer homogeneous polynomials, but sums of homogeneous 
polynomials^, as the following example shows^: 

Cil, 3) = ^ ■ (^1 + ^3) (70) 
Pi = 2 (s^ - 1) + l)s^-K ~s^ {s^ + l){s^ + 3s^-2)-E 
Ps ^ (6 - 1 + 11 s'^) ■ E^ + {s'^ ~ 1) (7 s^ + 12 s^ - 3) ■ K E^ 

+ {s^-l) {s^ + 3) (s4 + 2s2-l) {s^-l).EK^ 

+ {s^ - if {s^ - if . 

The two linear and cubic components Pi/3s^ and P^l3s^ are respectively solutions of 
the two linear differential operators: 

(3^4-7.^ + 14) . ii,4_9,2 + 4 



S (s2 + 1) (s2 _ 2) s2 (^2 _^ ly (g2 _ 2) (-1 + S^) 

-2 . + ^2 2 

(s2 - 1) s . iV " (s4 - 1)2 . TV ' 

■^^ + „4..4 .V. (71) 



s3(s4-l)2.iV ' s4(s4_i)3.7V 

N = 5^2 + 5 5I0 _^ 14 ,8 + 54 36 _^ 4g s4 _^ g2 _ 

^3 = 3s^'' + 15s^2 + 44si° + 98/ + 383/ + 415/ + 133s2-ll 
^2 = 19 s2" + 121 s^^ + 248 s^*^ - 408 s^^ - 974 s^^ + 2546 

+ 9597 s^ + 11440 + 6521 s"* + 1277 - 147 (72) 
= -27 s^"- 161 5^*^ + 240 3^6 + 5576 5^'' + 17854 s^^ + 28590 

tt This result can also be found in the Eqs. (3.22)-{3.35) of gl], which first show that C{N - \,N) 
is a homogeneous polynomial of E, K and the complete elliptic integral of the third kind EIi in the 
anisotropic case. Together with information from MontroU, Potts and Ward 1121 (note e.g. eq. A19) 
this means the same statement holds for C(Af — k,N) for k = 2,...,N. Reduction of Hi in the 
isotropic case then shows that C{N — k, N) is an inhomogeneous polynomial of E and K. 
^ Our results on the expressions of the C{N, M)'s are in agreement with those given , for N and M 
< 4, in 
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+ 30491 s^ + 19360 s^ + 8799 s^ + 1931 s^ - 333 

io 



An = -1792 5^"- 13136 s^^- 37568 s^^- 52256 s^''- 48848 



- 32576 s^° ~ 20720 s^ - 1568 s^ + 1600 s^ ~ 688 s^ + 192 

which are homomorphic to the first and third symmetric power of the hnear differential 
operator Le- 

equiv. SyirC'^LE), that is : Ly, ■ = W3 ■ Sym^{LE) 
Li equiv. Le, that is : Li ■ Qi = Wi ■ Le (73) 

where Q3 and W3 (resp. Qi and Wi) are linear differential operators of order three 
(rasp. one). The order six hnear differential operator corresponding to C(l, 3), that 
is the LCLM of Li and L3 is homomorphic to the LCLM of Le and Sym?{LE) ■ 

LiQLs equiv. LE®Sym^{LE) (74) 

Also note that for the horizontal, or vertical, correlations {N = or M ~ 0) 
one also has a homogeneous polynomials of E and K of degree zero. Let us consider 
for instance the simple correlation C(0, 1) : 



, Vl + S^ , (S-1)(S+1)V1 + S2 

C 0, n = 1/2^^ + 1/2-^ '-^ K 75 

s s 

The first term (of degree zero in E and K) is solution of an order one linear differential 
operator Zq, whereas the second term is solution of an order two linear differential 
operator li. 

lo ^ Ds + ,^ ^ (76) 



h = Ds 



2 



.S(l + s2)' 

{s^ -3)Ds 2 s*^ + 9 + 4 s2 + 1 



S(s2-1) ' (l + s2)2^2(g2_ 1X2- 



Up to a conjugation by (1 + s2)i/2^ ^^j^e order two linear differential operator li is an 
operator homomorphic to Le- 

(l + s2)-V2./^.(l + , 2)1/2 ^ (77) 



, (-4s2 + 3s4-3) s'^-s^ + ls^ 

Ds^ + ^- — — f • Ds 



(1 + S2)s(s2_l) ' (s2 _ 1)2 (1 + ^,2) g2 

with (1 + s2)-i/2 . . ^2)1/2 One actually finds that C(0, 1) is 

solution of the third order operator direct sum of Iq and h and is thus equivalent (up 
to conjugation by (1 + §2)^/2) ^j^g direct sum of la and Le- 

From the fact that the C{N, M)'s are actually always given by polynomial 
expressions sums of homogeneous polynomials in E and K, one easily deduces that 
the corresponding linear differential operators L^r m are homomorphic to direct sums 
of symmetric products of the second order linear differential operator l|4()|l. yielding 
generalizations of H36|l : 

Emm equiv. 0m 

5ym"(L£) (78) 

where for N — M odd, m is running as A^, + 1, iV + 2, • • • , M and for N — M even, 
as iV, A^ + 2, Af + 4, • • • , M, and where Sym"^{LE) = k) when m = 0. 

This structure is a consequence of the fact that the C{N, M)'s are given 
by polynomial expressions in E and K, instead of the rational or algebraic 
expressions in E and K which one could expect at first sight from the discrete 
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Painleve double recursions. This corresponds to quite remarkable identities and 
simplifications (factorizations, occurrence of perfect squares). From a less non-linear 
and more "Fuchsian" linear viewpoint, an explanation is the following. The non- 
diagonal C'{N, M) are determinants of holonomic functions, hence they are holonomic 
themselves. On the other hand, they are rational (or even algebraic expressions in E 
and K). Now, because the Galois group oi Le is SL{2, C), results from |^|^ show 
that expressions in E and K which are holonomic will have to be polynomial. 

Again one can check that all these linear differential operators L n m are Fuchsian 
differential operators with only three regular singular points t = 0, t — 1, t — oo. 
This is a straight consequence of the fact that these Latm's can be built as linear 
differential operators having polynomial solutions in E and K and thus, they inherited 
the three regular singular points t = 0, t = l, t — oo from the complete elliptic 
integrals E and K, and from the fact that the coefficients of the monomials E^ ■ 
are extremely simple rational expressions with no singularity except poles at s — Q 
(polynomial in s divided by powers of s). 

The results we got on the non-diagonal correlation functions C{N, M) are too 
numerous, and require too much space, to be given here (even if the final result is 
remarkably simple and elegant). However one sees the emergence of quite fascinating 
structures relating an infinite set of Fuchsian linear differential operators depending 
on two integers N and M (the Lnm's), with some quadratic double recursions that 
are nothing but discrete generalizations of Painleve, these structures being themselves 
closely linked with complete elliptic integrals. 



5. Backlund transformation and Malmquist Hamiltonian structure 

Let us recall that since the work of Malmquist it has been known that Painleve 
VI equation can be obtained from Hamilton equations 

dp dH dq dH 

with 

t{t-l)-H = qiq-l){q-t)p^ - Q{q) ■ p (80) 

+ (?i3 - ni) (713 - 712) (g - <) , where: Q(g) = 

(713 + rii) {q - 1) (g - <) + (rt3 - n^) q{q~t) ~ {ni + ^2) {q - 1) q 

With this structure, it follows that p is a rational function of t, q and qf. The 
Hamiltonian is the t— logarithmic derivative of the function r(t). The correlation 
functions C{N, N) being solutions of the sigma form of Painleve VI, one may find how 
the expressions of the two variables p and q (for which the Backlund transformations 
are birational) in the restricted case ni = N/2, 712 — (1 — N)/2, = (1 + N)/2 
and 714 = N/2 appear in terms of the elliptic integrals K and E. Considering the 
diagonal correlation function C(2, 2) taken as r(t) = t-^/^C(2, 2) one might expect, 
at first sight, to obtain the variables p and q as algebraic expressions in terms of E 
and K (and t). Remarkably, one obtains the surprising result that the two variables 
p and q are actually rational expressions*^ of E and K. For N = 2 one thus gets 

% Formulas expressing p and q as ratios of tau functions can be found in eqs. (5.42), (5.43) of 
Forrester and Witte 1361 . 
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two solutions, the simplest one being: 

{{t + i)E + (t-i)K) N^'' -n;,' 



P ■ 



2t {2E+{t-l)K) D^p'^ D^p^ 
t{2E+{t-l)K)-Nq 



{{t+l)E+{t-l)K) 



(1) 



(81) 



N^^^ = -{9t-l){t-lf ■ -2 {17 1 - 1) {t - 1) ■ EK 

+ {l + f - 34t) • E^ 

N^"^^ ^-{t-l)K^ -2EK + E'^ (82) 

= -3 {t -if -8{t- 1) EK + (-5 + t) E"^ 

^ (^t - if + 2 {t - if EK + {5t - 1) E'^ 

Nq ~ {3t ~ 11) {t - i f ■ K^ + 2 {t - 1) {3t^ ~ t + U) ■ EK 

+ {17 - 2t + 17) • E'^ 

One notes the homogeneous occurrence, in terms of degree, of E and K in these 
relations. The variables p and q have the rational parametrization of an algebraic 
curve. Obviously, the uniformization parameter similar to the one introduced in Sec. 
01 can be chosen as the ratio u = E/K (or E'/E) of the two elliptic integrals. 
One can then deduce that the parameter u is a solution of a Riccati differential 
equation. These results generalize straightforwardly to all the p, q associated with 
the C{N, Nys leading, remarkably, to rational functions of E and K and yielding 
rational parametrization for the corresponding algebraic curves between p and q. We 
have the same results in the variables a and a'. The expressions of the Backlund 
transformation corresponding to changing N into iV + 1 in terms of the variables p, 
q will be analyzed elsewhere. 



6. Conclusion 



The phenomenon of the existence of a one parameter family of solutions to Painleve VI 
equation has been presented in this paper by the study of the specific PVI equation 
which is satisfied by C{N, N) the diagonal correlation function of the Ising model. 
However the existence of such linear equations is a much larger phenomena and 
certainly holds for all PVI equations where the difference of any two of the parameters 
Vj is an integer because, in that case, there is a class of solutions which can be written 
as finite dimensional determinants whose elements are hypergeometric functions. 

Even though the existence of these Fuchsian differential equations follows from 
the general theorem on holonomic functions the specific form and properties of these 
equations is tedious to obtain. However, the expressions obtained for small N (via 
series computations) have been sufficient to guess the structure that is proved in 
sections 2 and 3. Moreover, using these initial computations, it has been possible to 
make a remarkably simple conjecture for the exponents which is in complete agreement 
with the local expansion of the Painleve VI equation at its singular points and this 
conjecture puts restrictions on the coefficients in the differential equations. 

In this paper we have obtained the Fuchsian equations by starting with the 
PVI equation. However the question can be reversed and we can ask what are the 
conditions on the Fuchsian equations which will lead to PVI equations. For second 
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order Fuchsian equations it would be sufficient to require that the exponents at the 
singularities agree with the exponents allowed by the local expansions of PVI. But for 
higher order equations the exponents do not fully specify the Fuchsian equation. The 
extra parameters which need to be specified are referred to as accessory parameters 
and only very specific accessory parameters will lead to Fuchsian solutions of PVI. 
The needed restrictions on these parameters are not known. 

The more general version of this is the question of determining whether or not a 
specific set of solutions to a Fuchsian equation will also satisfy some nonlinear equation 
(not necessarily PVI) . This is in some sense the original question asked by Jimbo and 
Miwa 0] and this is particularly important because, for C{N, TV), the nonlinear PVI is 
much simpler than the linear equations Lm.n- It was found in [ 1411151 [TBlll7| that the 
three and four particle contributions to the susceptibility of the Ising model, x, satisfy 
Fuchsian equations whose structure appears rather complicated and the question may 
be asked whether these functions, or their sum x, can also satisfy a nonlinear equation 
which might be simpler in appearance. 

Finally we remark that perhaps the most interesting discovery in this paper is 
that the operator Lmn are equivalent to the A^*'* symmetric power of the operator Le- 
This property extends to the operator Lh (which is isomorphic to Le)- One might 
wonder whether all solutions of the sigma form of Painleve VI that are also solutions 
of linear differential equations would be produced from symmetric powers oi Le by 
intertwinners. 
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